BAa4R-CONF-04/044 
SLAC-PUB- 10644 



Measurements of Branching Fractions and 
CP- Violating Asymmetries in S-Meson Decays to the 
Charmless Two-Body States K^tt+, K^K+, and 

The BABAR Collaboration 
August 17, 2004 

Abstract 

We present preliminary measurements of branching fractions and CP-violating 
asymmetries in decays of B mesons to two-body final states containing a K^. The 
results are based on a data sample of approximately 227 million T{4S) BE de- 
cays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at 
SLAC. We measure B{B+ K^tt+) = (26.0 ± 1.3 ± 1.0) x IQ-^, B{B+ K^K+) = 
(1.45l[]j^±0.11)xl0^6(< 2.35x10-6), and 5(5° ^ KO^O) ^ (l.l9t°;^°±0.13) x 10-^ 
where the first uncertainty is statistical and the second is systematic, and the upper 
limit is at a 90% confidence level. The significance of the B{B^ — > K^K^) and 
B{B^ — > K^K^) results are 3.5cr and 4.5cr, respectively, including systematic uncer- 
tainties. In addition, we obtain a measurement of the CP-violating asymmetry for the 
B"*" K^TT^ mode and we determine a 90% confidence-level interval for the asym- 
metry in the B+ K'^K+ mode: Acp{B+ ^ K'^tt+) = -0.087 ± 0.046 ± 0.010 and 
Acp{B+ K^K+) G [-0.43,0.68]. 
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1 Introduction 



The decays of B mesons into charmless hadronic final states provide important information 
for the study of CP violation. In particular, the study of the two-body decays B — > tttt, 
B Kn, and B — > KK provides crucial ingredients for measuring or constraining the values 
of the angles a and 7, defined as the following ratios of elements of the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix [1]: a = arg [-K.^K/fe/l/ ^1/* ] and 7 = arg [-F^^V^y l^^^V;;] . 
In this paper, we present measurements of the branching fractions of i?-meson decays to 
the charmless two-body final states K^7t+, K°K+, and K^K^. ^ For the B+ K°n+ and 
B^ — > K^K'^ modes we also report measurements of the direct CP asymmetries in the decay 
rates, 

_ T{B-^Klh-)-T{B+^Klh+) 

r {B- K^h-) + T{B+ ^ K^^h+) ' ^ ' 

where h = K,7r. 

Measurements of the rates and charge asymmetries for B Kn decays can be used 
to establish direct CP violation and to constrain the angle 7 [2]. The decay B^ K^ti^ 
is dominated by the h — > s-penguin process and in the Standard Model (SM) is expected 
to have Acp < 1% [3]. Thus, an observation of a sizable charge asymmetry could be an 
indication of non-SM contributions to the penguin-loop amplitude [3, 4]. The previously 
unobserved B KK decays proceed via penguin and M^-cxchange processes similar to 
those in B^ — > tt^tt^ and can help in the determination of a in the measurement of time- 
dependent CP asymmetries in B^ n'^n' [5]. Measurements of the branching fractions for 
these decay modes also provide important information regarding rescattering processes [6] . 



2 The BaBAR detector and dataset 

The measurements presented in this paper are based on data collected with the BABAR 
detector [7] at the PEP-II asymmetric-energy e'^e~ collider [8], located at the Stanford Linear 
Accelerator Center. The sample consists of 226.6 ± 2.5 miUion BB pairs produced at the 
T{4:S) resonance ( "on- resonance" ) , which corresponds to an integrated luminosity of about 
205 fb""^. An additional 16fb~^ of data recorded at an e'^e~ center-of-mass (CM) energy 
approximately 40MeV below the T{4S) resonance ("off-resonance") is used for background 
studies. 

The BABAR detector is described in detail in [7]. Charged-particle momenta are measured 
in a tracking system consisting of a five-layer, double-sided silicon vertex detector and a 40- 
layer drift chamber (DCH), which operate in a solenoidal magnetic field of 1.5 T. Particles 
are identified as pions or kaons based on the Cherenkov angle measured with a detector 
of internally refiected Cherenkov light (DIRC). The direction and energy of photons are 
determined from the energy deposits in a segmented CsI(Tl) electromagnetic calorimeter 
(EMC). 

^Unless explicitly stated otherwise, charge-conjugate decay modes are assumed throughout this paper 
and branching fractions are averaged accordingly. 
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3 Analysis method 



Hadronic events are selected on the basis of charged-particle multiplicity and event topology. 
We reconstruct i?-meson candidates decaying to K^X, where X can be tt"*", K~ or K^. The 
candidates are reconstructed in the mode ^ ^ n'^n'. 

The following selection criteria are applied to the candidate S-decay products. Charged 
tracks are required to be within the tracking fiducial volume and to have at least 12 DCH 
hits and a minimum transverse momentum of O.lGeV/c. Tracks that are not iT^-decay 
products are required to originate from the interaction point, to be associated with at least 
six Cherenkov photons in the DIRC, and to have a Cherenkov angle within 4 cr of the expected 
value for a pion or kaon hypothesis. Candidate mesons are reconstructed from pairs of 
oppositely charged tracks that are consistent with originating from a common vertex, have 
an invariant mass within ±11.2 McV/c^ of the nominal mass, and have a measured proper 
decay time greater than five times its uncertainty. 

The S-meson candidate is characterized by two nearly uncorrelated kinematic variables: 

the energy-substituted mass, uies = \J {s/2 + pj ■ p^)'^ / Ef — p\, and the energy difference, 
A£' = E*Q — y/s/2, where the subscripts i and B refer to the initial e^e~ system and the B 
candidate, respectively. The asterisk denotes the T{AS) rest frame (CM frame), and ^/s is 
the total CM energy. The pion mass is assigned to all charged particles in the calculation 
of E^. For 5° K^K^ candidates, we require |AE| < 0.1 GcV, while for B+ K°h+ 
candidates {h = tt, K) we require —0.115 < AE < 0.075 GeV. The interval is asymmetric in 
order to select both 5+ — > K^n'^ and 5+ — > K^K'^ decays with nearly 100% efficiency. The 
AE distribution peaks near zero in the modes not containing charged kaons. In S+ — > K'^K'^ 
decays, the AE distribution peaks at —45 MeV as a result of the assignment of the pion mass 
to the charged kaon candidate. The distribution of itles peaks at the B mass in all modes. We 
impose a loose selection of 5.20 < tties < 5.29 GeV/c^ that includes a background-dominated 
region used to estimate the level of backgrounds in the signal region. 

Simulated events [9], off- resonance data, and events in on- resonance mss- and AE- 
sideband regions arc used to study backgrounds. The contribution from other S-meson 
decays is found to be negligible. The primary source of background are random combina- 
tions of tracks and neutral clusters produced in e+e~ — > qq events, where q is a u, d, s, or c 
quark. In the CM frame, these "combinatorial" events are characterized by a jet-like struc- 
ture, in contrast to the more uniformly distributed decays of B mesons produced in T{4S) 
decays. We suppress combinatorial backgrounds by exploiting this topological difference 
through a selection based on 6g, the angle between the sphericity axis of the 5-candidate 
decay products and the sphericity axis of the remaining particles in the event. In the CM 
frame, this selection is | cos 6*5! < 0.8. To discriminate further between signal decays and 
combinatorial backgrounds, we employ a Fisher discriminant, J^. We define as an op- 
timized linear combination of ^jP* and ^iP* cos^ 9* [10], where p* is the momentum of 
particle i and 6* is the angle between its momentum and the thrust axis of the 5-candidate 
decay products, both calculated in the CM frame. The sums arc over all particles in the 
event except for the i?-candidate decay products. The difference between signal and back- 
ground distributions of is present in the probability density functions (PDF's) that model 
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T in the fits to the data sample, which we describe below. 

Signal yields and charge asymmetries are determined from unbinned extended maximum- 
likelihood fits. The extended likehhood for a sample of N K'^X candidates is defined as 



N 



exp -^n^m 



(2) 



where Vi{xj; oti) is the PDF of a signal or background category i, evaluated at the values of the 
measured variables Xj of candidate j. The sum is over the set of categories. The parameters 
di determine the expected distributions of measured variables in each category, and rii are 
the yields being determined in the fit. The probability coefficients Ni are defined separately 
for each mode. We perform separate fits to the two samples of B candidates: B^ — > K^K^ 
and B+ K^h+ {h ^n,K). In the fit to the "neutral-5" (5° ^ K^K^) sample we include 
two categories, signal and background. The probability coefficients Ni are set equal to the 
yields (iVj = rii); the yield in each category is obtained by maximizing the likelihood. In the 
fit to the "charged-i?" (5+ — > K^h^) sample, we fit simultaneously two signal categories, 
B^ — > ir°7r+ and B^ — > K^K^, and two corresponding background categories. In addition, 
the probability coefficient for each category i is given by A^j = rij (1 — qjAi)i where rij is the 
total yield, summed over charge states, Ai is the charge asymmetry, and Qj is the charge of 
the B candidate. The total yields and charge asymmetries are determined by maximizing 
C. 

As the input variables in the fit are nearly uncorrelated, the PDF in the likelihood 
function, Vi{xj;di), is constructed as the product of the individual PDF's of the input 
variables Xj. In both fits, the set of input variables contains mEs, AE", and J^. In the 
fit to the B^ K^h^ sample, we include also the normalized Cherenkov-angle residuals 
{Oc — 0^) jue^ and (6'c — 6'^) l^e^-, where Qc is the measured Cherenkov angle of the primary 
daughter oq^ is its measurement uncertainty, and Q1 (9^) is the expected Cherenkov 
angle for a pion (kaon) hypothesis. The quantities ae^, 9'^, and 9^ are measured separately 
for negatively and positively charged pions and kaons from a control sample of D° — > K'~7r'^ 
decays originating from D*^ decays. 

The parameterization of the PDF's is determined from data and Monte Carlo-simulated 
(MC) samples. Some PDF parameters are free to vary in the fit as explained below. The 
signal niES PDF's for S+ — > K^h'^ and B° — > K^K'^ are derived from signal MC samples. 
The shape of the distribution is modeled as a Gaussian function with an asymmetric variance 
and a low-side power tail in the K'^h^ fit, while in the K^K^ fit it is parameterized as a linear 
combination of two Gaussian functions ("double- Gaussian" function). The mean value of the 
signal mES distribution is a free parameter in the S+ — > fit, as the sample of candidates 

is sufficiently large. To describe the background mES PDF, we use an empirical threshold 
function [11]. The shape parameter of this function is a free parameter in the B^ K^h^ 
fit, while in the B^ K^K^ fit it is determined from the on-resonance A£J-sideband region, 
defined by the selection 0.1 < \ AE\ < 0.3 GeV. 

The signal AE PDF's are derived from MC samples and modeled as double-Gaussian 
functions for both modes. The signal AE distribution is expected to be centered near zero 



11 



Tabic 1: Summary of results for numbers of seleeted K^X candidates A^, total detection effi- 
ciencies e, fitted signal yields Ns, statistical significances S including systematic uncertainties, 
charge-averaged branching fractions B, and charge asymmetries Acp- The efficiencies include the 
branching fractions of intermediate states (K^ Kg tt+tt"). Branching fractions are calculated 
assuming equal rates for the T(4S') — > B^B^ and T{4S) B^B~ processes. The upper limit 
for the K^K^ branching fraction corresponds to a 90% confidence-level (C.L.), and the central 
value is given in parentheses. The 90% C.L. asymmetry interval is given for K^K'^, including the 
systematic uncertainty. The 90% C.L. asymmetry interval for K'^tt'^ is [—0.164,-0.010]. 



Mode 


N 


s{%) 




Ns 




S{a) 


B(10-«) 


Acp 


K°Tr+ 
K°K+ 


20441 


12.6 ±0.3 
12.5 ±0.3 


744 
41 


+37 
-36 

+15 

-13 


+21 
-17 

+3 
-2 


3.5 


26.0 ± 1.3 ± 1.0 
< 2.35 (1.45tg j^ ± 0.11) 


-0.087 ±0.046 ±0.010 
[-0.43,0.68] 




I!):-')!) 


(S.T) ± 0.5 


2:-!.() 


+ 7.7 

-(i.T 


-1.9 
-2.0 


4.r, 







for the K^n^ mode, while for K^K^ candidates, the mean of AE' is shifted because the pion 
mass is assumed for the charged track. The shift in is 



(A^) = -7boost X 




where p is the momentum of the track, and and Mk are the nominal values of the 
pion mass and the kaon mass, respectively. The Ai? mean value is a free parameter in 
the — > K^h^ fit, while in the K^K^ fit it is determined through a comparison 

of the values obtained in the two MC samples with the value obtained in the fit to the 
— > K^h^ data sample. The background A£' distribution is modeled as a second- and 
first-degree polynomial function for the charged-S and neutral-5 modes, respectively. The 
polynomial coefficients are determined from on-resonance events in the mEs-sideband region, 
defined by the selection 5.20 < < 5.26GeV/c^. 

In both modes, the signal T distribution is modeled as a Gaussian function with an 
asymmetric variance [12]. Its parameters are free to vary in the K^h^ fit; in the K^K^ fit, 
they are determined from the MC sample. The background distribution is parameterized 
as a double-Gaussian function with its parameters left free to vary in both fits. 

In the charged- i? modes, the normalized Cherenkov-angle residuals are modeled as double- 
Gaussian functions; the PDF parameters are taken from the D* control sample and they are 
determined separately for tt"*", tt", , and K~ tracks as a function of momentum and polar 
angle. 

4 Physics results and systematic uncertainties 

The results of the two maximum-likelihood fits are summarized in Table 1. The K^K^ final 
state is an equal admixture of K^gK^ and K^K^. We therefore use a 50% probability for the 
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Figure 1: Distributions of mES and /\E for (a,b) B+ i^^7r+, (c,d) B+ K^K+ and (e,f) 
— > KgKg (histograms) after background subtraction (see text). Projections of the fit PDF's 
are overlaid (sohd curves). 13 
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Figure 2: Distributions of iriES and AE for (a,b) K^h^ and (c,d) K^K^ (histograms) 

after signal subtraction (see text). Projections of the fit PDF's are overlaid (solid curves). 



14 



K^K^ to decay as K^K^ in computing the — > K^K^ branching fraction. We also use the 
current world averages for B{Kg — > tt+tt") [13]. 

Figure 1 shows background-subtracted distributions of tties and AE for — > K^h^ 
and B^ K^K^ candidates. The background subtraction is performed by weighting events 
using the sVlot technique described in Ref . [14] . The shape of the resulting distribution can 
be compared with the PDF used in the full fit. We find good agreement in both variables 
for K^h^ and K^K^ candidates. The corresponding signal-subtracted distributions of mss 
and AE are shown in Figure 2. 

The signal significance is defined as the square root of the difference in —2 In C between 
the best fit and the null-signal hypothesis. The upper limit on the signal yield for a given 

mode i is defined as the value of n\^^ for which J^"' Cia^drii/ C-^aa^drii = 0.9, where Cma^ 
is the likelihood as a function of n,, maximized with respect to the remaining fit parameters. 
Branching-fraction upper limits are then calculated by increasing the signal-yield upper limit 
and reducing the efficiency by their respective systematic uncertainties. 

Systematic uncertainties in the signal yields arise primarily from imperfect knowledge of 
the PDF parameterizations. Such systematic errors are evaluated either by varying the PDF 
parameters by their measured (Icr) uncertainties or by substituting alternate parameteri- 
zations. The significance of each signal yield with the systematic uncertainties included is 
evaluated by imposing simultaneously in the fit all of the systematic variations of the PDF's 
that lower that signal yield. Systematic uncertainties in the efficiency include uncertainties 
in tracking and Kg reconstruction. 

In the K^h^ sample, the dominant systematic uncertainty is that associated with the 
signal mEs shape, leading to a systematic error of Ij^g (lii?) events, and that associated with 
the signal AE resolution, leading to a systematic error of (to 7) events in the B^ K'^tt'^ 
[K^K^) mode. The significance of the B^ K^K^ observation with (without) the sys- 
tematic uncertainty included is 3.5 a (3.7 cr). In the B^ — > K^K^ mode, the main systematic 
uncertainty is that associated with the background toes distribution (±0.8 events). There 
is also a sizable systematic contribution from the imperfect agreement between MC and 
data samples (±1.4 events) and, in the branching-fraction extraction, the uncertainty in the 
K^g efficiency (6.0%). The significance of the B^ K^K^ observation with (without) the 
systematic uncertainty included is 4.5(7 (4.8 cr). The systematic uncertainties in the charge 
asymmetries in the B^ — > K'^n'^ mode are evaluated by adding in quadrature the contri- 
butions from PDF variations and the upper limit on intrinsic charge bias in the detector 
(±0.010). 

5 Summary 

We present preliminary measurements of the branching fraction and the CP- violating charge 
asymmetry in the B^ — > K^t:'^ decay, and a preliminary measurement of the branching 
fraction of the S° — > K^K^ decay. No evidence of direct CP violation in the B^ — > K^'k'^ 
mode is observed. The B^ K^K^ measurement constitutes the first observation of this 
decay channel: the probability of obtaining our result assuming the null-signal hypothesis 
is 3.4 X 10~®. We have also searched for the B^ — > K^K~^ decay and set an upper limit on 
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its branching fraction at 2.35 x 10~^ at the 90% C.L. The branching-fraction measurements 
reported here are consistent with previous measurements of the same quantities [15, 17, 18, 
16], but are extracted from a data sample larger by a factor of 2.6. 
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